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^ C"j Violation of Lorentz invariance and CPT symmetry is a predicted phenomenon of 

Planck-scale physics. Various types of data are analyzed to search for Lorentz violation 
under the Standard-Model Extension (SME) framework, including neutrino oscillation 
data. MiniBooNE is a short-baseline neutrino oscillation experiment at Fermilab. The 
measured excesses from MiniBooNE cannot be reconciled within the neutrino Standard 
Model (fSM); thus it might be a signal of new physics, such as Lorentz violation. We 
have analyzed the sidereal time dependence of MiniBooNE data for signals of the possible 
breakdown of Lorentz invariance in neutrinos. In this brief review, we introduce Lorentz 
violation, the neutrino sector of the SME, and the analysis of short- baseline neutrino 
1 oscillation experiments. We then present the results of the search for Lorentz violation 

f****) in MiniBooNE data. This review is based on the published result . 
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1. Spontaneous Lorentz symmetry breaking (SLSB) 

Every fundamental symmetry needs to be tested, including Lorentz symmetry. The 
breakdown of Lorentz invariance naturally arises in different scenarios of physics at 
the Planck scale. For this reason the expected scale of Lorentz- violating phenomena 
is more than the Planck mass mp ~ 10 19 GeV, or in other words, Lorentz viola- 
tion is expected to be suppressed until at least ~ 10~ 19 GeV in our energy scale. 
Lorentz symmetry is a fundamental symmetry both in quantum field theory and 
general relativity; the consequence of its violation would be enormous, and it seems 
it is impossible to establish a self-consistent theory with Lorentz violation. However, 
Lorentz violation can be incorporated into existing theories by spontaneous break- 
ing. In this way, Lorentz-violating terms in the Lagrangian do not conflict with the 
Standard Model (SM). 

There are a number of models for spontaneous Lorentz symmetry breaking 
(SLSB), 2 although the basic idea is the same for all. In the SM, the mass terms 
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arise from the spontaneous symmetry breaking (SSB) triggered by a Higgs field 
with a nonzero vacuum expectation value. In this way, mass terms are dynamically 
generated and they do not exist before SSB. Consequently, the SM does not have 
to violate gauge symmetry where mass terms would necessarily do so. 

Figure 1 illustrates this situation. 3 The theory starts from perfect symmetry and 
the vacuum is true null space (Fig. la). After the SSB, the vacuum is saturated by a 
scalar Higgs field <f>, and the particles obtain mass terms (Fig. lb). This idea can be 
extended to a vector field (Fig. lc). The ultra- high-energy theories, such as Planck- 
scale physics theories, have many Lorentz vector fields (or, more generally, Lorentz 
tensor fields). When the universe cools, if any of them acquire nonzero vacuum 
expectation values, then the vacuum can be saturated with vector fields. (Fig. Id). 
Note that, theoretically, SLSB is conceived to occur earlier than the SSB of the SM, 
unlike this cartoon. Such vector fields are the background fields of the universe, and 
they are fixed in space; couplings with the SM fields generate interaction terms in 
the vacuum. 

L = ixpj^tp + rmpip + ipj^ip + ip^c?" d„ip H . 

In this expression, the coefficients a M and c^ v represent vacuum expectation 
values of vector and tensor fields, and they correspond to background fields that fill 
the universe. The crucial observation is that, since these Lorentz tensors are fixed 
in space and time, they cause direction- dependent physics. In particular, rotation 
of the Earth (period 86164.1 sec) causes sidereal time dependent physics for any 
terrestrial measurement if the SM fields couple with Lorentz- violating background 
fields. Therefore, the smoking gun of Lorentz violation is to find a sidereal time 
dependence in any physics observable. 

2. What is Lorentz and CPT violation? 

We introduce Lorentz violation as coupling terms between ordinary SM fields and 
background fields in the universe. They are Lorentz scalars of the coordinate trans- 
formation; however since background fields are fixed in space, motion of the SM 
particles generates coordinate-dependent physics. 

The situation is illustrated in Figure 2. The top cartoon (Fig. 2a) shows our 
setting: a SM particle is moving in two-dimensional space, from bottom to top, as 
seen by the local observer (Einstein). The space is filled with a hypothetical Lorentz- 
violating background field, (depicted by arrows). There are two ways to move 
this particle from left to right for the local observer: Particle Lorentz transformation 
and Observer Lorentz transformation. 

2.1. Particle Lorentz transformation 

The first one is Particle Lorentz transformation, where the motion of a SM particle 
is actively transformed in the fixed coordinate system (Fig. 2b). Since the back- 
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Fig. 1. An illustration of spontaneous symmetry breaking (SSB). 

ground field is unchanged, as a consequence, a coupling between the SM particle 
and the background field is not preserved; therefore, one can see Lorentz violation. 
In other words, Lorentz violation generally means the Particle Lorentz violation, 
and it implies direction-dependent physics of SM particles in the fixed coordinate 
space. 



2.2. Observer Lorentz transformation 

The second way is Observer Lorentz transformation, where the coordinate is in- 
versely transformed (Fig. 2c). In this cartoon, if Einstein (the local observer) turns 
his neck 90° counterclockwise, the SM particle moves from left to right for the 
local observer. However, the background field is also transformed to the new coor- 
dinate system leaving the coupling with the SM particle unchanged. One cannot 
generate Lorentz violation by Observer Lorentz transformation, because Observer 
Lorentz transformation only corresponds to a coordinate transformation. In other 
words, coordinate transformations preserve the Lorentz-violating effect, and every 
observers agree with the same Lorentz-violating effect. 
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b) Particle Lorentz 
transformation 
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Lorentz violation! 



c) Observer Lorentz 
transformation 




Fig. 2. An illustration of Particle Lorentz transformation and Observer Lorentz transformation. 



2.3. CPT violation 

There is a close link between Lorentz symmetry and CPT symmetry. Here, "CPT" 
represents the combination of charge transformation (C), parity transformation (P), 
and time reversal (T). It is known that none of these, taken individually, is a sym- 
metry of the SM, but from the CPT theorem 4 we expect that their combination 
is a perfect symmetry. Since the CPT theorem is based on Lorentz symmetry, one 
can expect CPT violation when Lorentz invariance is broken. This is manifest in 
the appearance of CPT-odd terms in the Lagrangian, which appear as a subset of 
the terms that break Lorentz invariance. The phase of CPT transformation is re- 
lated to the number of Lorentz indices, i, transforming under the Particle Lorentz 
transformation. 



CPT phase = (-1)' . 

In the SM any Lagrangian is the linear sum of CPT-evcn, or phasc= +1 terms. 
This is the main consequence of the CPT theorem, and this is why any Lagrangian 
is CPT invariant. However, if the theory includes Lorentz violation, it is possible 
that i can be odd. When i is odd number, Lorentz violation causes CPT violation, 
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and this is called "CPT-odd" Lorentz violation. On the other hand, if i is an even 
number, the phase of CPT transformation is even, and the theory does not violate 
CPT, even though it contains Lorentz violation. Such Lorentz violation is called 
"CPT-even" Lorentz violation. CPT-odd and CPT-even Lorentz violation differ 
only in the number of their Lorentz indices. 

• coefficients of CPT-odd Lorentz- violating terms (a M , g^ vX , • • • ) 

• coefficients of CPT-even Lorentz- violating terms {c^ v ', K MiyAcr , • • • ) 

Note that the interactive quantum field theory necessarily violates Lorentz sym- 
metry if CPT symmetry is not preserved. 5 This general theory is consistent with 
the argument here and Standard-Model Extension (SME), which we discuss in the 
next section. 

3. Analysis of Lorentz violation 

Lorentz violation is realized as a coupling of SM particles and background fields. 
To specify the components of the Lorentz- violating vector or tensor fields, the coor- 
dinate system must be specified. Then the general Lagrangian, including all possi- 
ble Lorentz-violating terms, is prepared. Finally, using this Lagrangian, observable 
physical quantities can be identified. 

3.1. The Sun-centered coordinate system 

The choice of the coordinate system is arbitrary, since Particle Lorentz violation in 
one coordinate system is preserved in another coordinate system through Observer 
Lorentz transformation. Nevertheless, in order to compare experimental results from 
different experiments in a physically meaningful way, a common frame should be 
used. For this purpose, we need a universal coordinate system (Fig. 3). 6 The uni- 
versal coordinate system is required to be reasonably inertial in our timescale. The 
Sun-centered coordinated system is just such a coordinate system (Fig. 3a). Here, 
the rotation axis of the Earth aligns with the orbital axis by tilting the orbital plane 
by 23.4°, defining the Z-axis. The X-axis points towards the vernal equinox, and 
the Y-axis completes the right-handed triad. Obviously, we assume Lorentz-violating 
field is uniform at least the scale of the solar system. This can be justified in many 
ways, for example, we know the weak and the electromagnetic laws are same in far 
stars through the observation, indicating Lorentz-violating fields are also uniform 
in these scales if they are arisen through the spontaneous breaking process. Note, 
Sun-centered is more suitable than galaxy-centered, because although the galac- 
tic rotation is faster, the galactic rotation takes too long for human observation to 
change the direction, and consequently it cannot help test the violation of rotational 
symmetry. Then the location of the experiment is specified by the Earth-centered 
coordinate system (Fig. 3b). Here, the x-axis points south, the y-axis points east, 
and the z-axis points to the sky from the site of the experiment. Finally, local polar 
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Fig. 3. The coordinate system used by this analysis; (a) first, the motion of the Earth is described 
in Sun-centered coordinates, (b) then the local coordinates of the experiment site are described in 
Earth-centered coordinates, (c) finally, the direction of the neutrino beam is described in the local 
polar coordinate system, (d) The time zero is defined when the experiment site is at midnight near 
the autumnal equinox, in other words, when the large "Y" and small "y" axes almost align. 

coordinates specify the direction of the beam (Fig. 3c). The time zero of the side- 
real time is defined as being the position of the experiment at midnight near the 
autumnal equinox (Fig. 3d). 

3.2. Standard-Model Extension (SME) 

For the general search for Lorentz violation, the Standard-Model Extension 
(SME) 7:8 is widely used by the community. Various types of data are analyzed 
under SME, 10 ' 11 including neutrino oscillation data. 6 ' 12 ' 13 ' 14 ' 15 The SME is a self- 
consistent effective field theory including Particle Lorentz violation. In principle, 
SME is an infinite series of all types of interactions, but many analyses are lim- 
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ited to the renormalizable sector, called the minimal SME. For example, under the 
minimal SME, the effective Lagrangian for neutrinos can be written as, 9 



1 - <-> 

C = -itpAT^ D M ip B - ipAM AB i: B + h.c, (1) 

T AB = Y&AB + CABlv + d ABl^H + e AB + «/ab75 + ^1b^, (2) 

Mab = m, A B + im 5 ABj5 + o-abIh + b AB + 2 Hl AB a nv ( 3 ) 

Here, the AB subscripts represent Majorana basis flavor space (6x6 for con- 
vention). The first term of Eq. 2 and the first and second terms of Eq. 3 are the 
only nonzero terms in the SM, and the rest of the terms are from the SME. As we 
see, these SME coefficients can be classified into two groups (Sec. 2.3), namely e t AB , 
Jabi 9ab > a ABi an d b AB which are CPT-odd SME coefficients, and d% B , d A v B , and 
H 4 „ which are CPT-even SME coefficients. 



3.3. Lorentz-violating neutrino oscillations 

Once we have a suitable formalism, such as the SME, we are ready to write down 
physical observables. The effective Hamiltonian relevant for the v — v oscillations 
can be written, 9 



(h cS )ab ~ rp\{a L yPii - {c L y v PnP»]ab (4) 

Here, (ol)^ = (a + b)^ b and (cl)^ = (c + d)^. Note that we drop the neutrino 
mass term since the standard neutrino mass term is negligible for short baseline 
neutrino oscillation experiments, such as MiniBooNE. 

Solutions of this Hamiltonian have very rich physics, but for our purpose, we 
restrict ourselves to short- baseline — v e (v^ — v e ) oscillation phenomena. By 
assuming the baseline is short enough compared to the oscillation length, the — u e 
oscillation probability can be written as follows, 16 



P v ^ Ve ^ I ( C W + (A)e M sinu;©T© + (A c )en cosw©T© 

+ (£ s )e AI sin2w©T© + (B c ) ell cos 2w©T© | 2 . (5) 

Here, w© stands for the sidereal time angular frequency (w© = m ^ 4 1 rad/s), 
as opposed to the solar time angular frequency (ojq — 86 4q rad/s). The neutrino 
oscillation is described by the function of the sidereal time T© , with the sidereal time 
independent amplitude (C) eM , and the sidereal time dependent amplitudes, (*4 s ) eAt , 
(A c ) ei j,, (B s ) ef _ l , and (B c ) ell . Therefore, an analysis of Lorentz and CPT violation in 
neutrino oscillation data involves fitting the data with Eq. 5 to find nonzero sidereal 
time dependent amplitudes. 
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In terms of the coefficients for Lorentz violation, these amplitudes are explicitly 
given by, 



(B s ) e , 



(C^)^ + E(C^) e , 
(Af\ ti + E(A^) 



efj, 



(6) 
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N*N*)((c L )t* (c L ) Y J) - 2N x N Y (c L ) XY 



Here, the N x , A r , and A z are the direction vectors of the neutrino beam 
in the Sun-centered coordinates (Sec. 3.1). The components are described with a 
co- latitude \ of detector location in the Earth-centered system (Fig. 3b), and the 
zenith and azimuthal angles 9 and <f> of the local beam system. (Fig. 3c): 









■( 







cos x sin 8 cos 4> + sin \ cos 9 

sin 9 sin <j> 
- sin x sin 9 cos </> + cos x cos # , 



(8) 



For the antineutrino oscillation analysis, one needs to switch the sign of 
according to CPT-odd nature of CPT-odd coefficients (a^ — > ~a L ). 

In the reality of the analysis, fitting five parameters using Eq. 5 is not easy. 
Therefore, we also consider the following three-parameter model, by setting (B s ) efl 
and (B c ) efl to be zero by hand. This model, Eq. 9, can be motivated, for example, 
by assuming nature only has CPT-odd SME coefficients. 



(hey 



\(C) eft + (A s )e^ SmW e T e + (A c ) efi cosw ffi T ( 



(9) 
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Fig. 4. L-E diagram with the uSM (two straight dotted lines) and the Puma model (two dashed 
and solid curves). 



3.4. Lorentz violation as an alternative neutrino oscillation model 

Because of the unconventional energy dependence of Lorentz-violating terms in the 
Hamiltonian (E° and E 1 ), naively, its energy dependence on neutrino oscillations 
is different from one expected from the three massive neutrino model (so-called 
i^SM). However, it is also possible to "mimic" neutrino mass- like energy dependence 
(E^ 1 ) using Lorentz violating terms on ly. 17 < 18 ' 19 ' 20 There is a chance that such 
types of models could be correct, because we currently have some tensions in the 
world neutrino oscillation data. For this purpose, it would be helpful to show the 
phase space of neutrino oscillations in a model-independent way. The L-E diagram 
(Fig. 4) shows world's neutrino oscillation experiments mapped with their energy 
and baseline. 19 

The curves in Figure 4 represent the oscillation length. For example, massive 
neutrino oscillation solutions (=L/E oscillatory dependence) are represented by 
the line L oc E. Here, data are consistent with two L/E neutrino oscillations, the 
v e disappearance measurement at the KamLAND experiment (2 to 8 MeV), and 
the U/j, and disappearance measurements at the long-baseline and atmospheric 
neutrino experiments (0.3 to 10 GeV). Therefore, we know there are at least two 
segments with L oc E on the L-E diagram. Nevertheless, our knowledge outside of 
these segments is limited. There are proposed models, such as the Puma model, 19 ' 20 
which have L/E oscillatory dependencies in these energy ranges. So, here, the mod- 



August 14, 2012 0:36 WSPC/INSTRUCTION FILE LVreviewV6 



10 T. Katori 




primary beam secondary beam tertiary beam 
► ► » 

(protons) (mesons) (neutrinos) 



Fig. 5. An overview of the MiniBooNE experiment. The top left picture shows the Fermilab site, 
including the Booster. The top right picture is the tank of the MiniBooNE detector. The bottom 
cartoon shows the sketch of the Booster Neutrino Beamline (BNB). 

els are consistent with current data, a but outside of these energy ranges they have 
completely different dependencies. These alternative models are interesting because 
they have a chance to reproduce short-baseline anomalies, such as the MiniBooNE 
oscillation signals, 24,25 which we discuss in the next sections. 

4. The MiniBooNE experiment 

The MiniBooNE experiment is a short-baseline neutrino oscillation experiment at 
Fermilab, USA (2002-2012). Its primary goal is to find — > v e (v^ — > v e ) oscil- 
lations with an ~800 (600) MeV neutrino (antineutrino) beam with an ~500 m 
baseline. Figure 5 shows the overview of the MiniBooNE experiment. 26 

4.1. Booster neutrino beamline (BNB) 

MiniBooNE uses neutrinos (antineutrinos) from the Booster neutrino beamline 
(BNB), 27 which is illustrated in Figure 5, bottom. The 8 GeV protons, the "pri- 

a Recent reactor neutrino disappearance oscillation results do not support Puma model 21 > 22 > 23 
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mary" beam, are extracted from the Booster and steered to collide with the beryl- 
lium target in the magnetic focusing horn. The collision of protons and the target 
makes a shower of mesons, the "secondary" beam; and the toroidal field created by 
the horn focuses tt + (tt~) for neutrino (antineutrino) mode. At the same time it 
defocuses ir~ (ir + ), which create backgrounds. The decay-in-flight of tt + (ir~) create 
(vp), the "tertiary beam". The consequent muon neutrinos (muon antineutrinos) 
are a wideband beams peaked in around 800 (600) MeV. 

4.2. MiniBooNE detector 

The MiniBooNE detector is located 541 m to the north of the target. 28 The detec- 
tor is a 12.2 m diameter spherical Cherenkov detector, filled with mineral oil. An 
optically separated inner black region is covered with 1,280 8- inch PMTs, and an 
outer white region has 240 8-inch PMTs which act as a veto (Fig. 5, top right). The 
black inner cover helps to reduce reflections, so that one can reconstruct particles 
from the Cherenkov light more precisely; the outer white cover helps to enhance 
reflections, so that a smaller number of veto PMTs can cover a larger area. 

4.3. Events in detector 

The time and charge information of the Cherenkov ring from the charged particle is 
used to estimate particle type, energy, and direction. 29 For example, an electron-like 
track is characterized by a fuzzy-edged Cherenkov ring, compared with a sharp- 
edged muon-like Cherenkov ring. Based on particle type hypothesis, the track fit- 
ter estimates a particle energy and direction. Figure 6 shows typical particles and 
their characteristic tracks, Cherenkov rings, and event candidates from the event 
display. 26 

A variety of track fitters are developed to measure the kinematics of specific types 
of interactions. 30 > 31 > 32 > 33 Among them, the most important reaction for the oscilla- 
tion analysis is the charged current quasielastic (CCQE) interaction, 34 ' 35 which is 
characterized by one outgoing charged lepton (at the BNB energy, protons seldom 
exceed a Cherenkov threshold of ~350 MeV kinetic energy). If a charged lepton 
is detected from the CCQE interaction, one can reconstruct the neutrino energy, 
E® E , by assuming the target nucleon is at rest and the interaction type is truly 



Here, M n , M p , and m M are the neutron, proton, and muon masses, is the 
total muon energy, 9^ is the muon scattering angle, and B is the binding energy of 
carbon. Ability to reconstruct neutrino energy is essential for neutrino oscillation 
physics, since neutrino oscillations are function of neutrino energy. 



— V ' o < 

CCQE 34 (QE assumption). 



E# E = 



2(M W - B)E ll - ((M n - Bf + ml- M 2 ) 



(10) 
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Interaction track Cherenkov Candidate 




Fig. 6. The particle type and characteristics. From left to right, interaction types, characteristics 
of tracks and Cherenkov rings, and event displays of candidate events. 



It is interesting to note that many of the neutrino interaction cross sections 
measured by MiniBooNE are at a higher rate and harder spectrum than historically 
known values and disagree with interaction models tuned with electron scattering 
data. 36 This fact triggered the development of a new class of neutrino interaction 
models, 37,38 ' 39,40 ' 41, 42,43 mostly by including nucleon correlations. These new models 
even question how to reconstruct neutrino energy 44 ' 45 with the QE assumption, 
traditionally done in all Cherenkov-type detectors. Therefore, similar to other fields 
(e.g., cosmology), the further study of neutrino physics just increases the number 
of mysteries! 



4.4. Oscillation analysis 

The signature of the v u — ¥ v e (y^ — ¥ v e ) oscillation is the single, isolated electron-like 
Cherenkov ring produced by the CCQE interaction. 



oscillation _ 

i/fj, — > v e + n^ e +p , 

_ oscillation _ _i_ , 

v u — > v e + p ->• e~ + n . 



Note, since MiniBooNE was not magnetized, electrons and positrons were not 
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distinguished. Thus, the analysis of neutrino mode and antineutrino mode is rea- 
sonably parallel, except for some differences in handling backgrounds. 

There are two backgrounds which contribute equally to our signals. The first 
class is the "misID" , and this is dominated by a single gamma ray from the neutral 
current channels, such as radiative A decay and n° production, where only one 
gamma ray is detected. It is essential to constrain our misID background predictions 
on those channels using measurements of controlled samples. For this purpose, we 
measured the neutral current 7r° production rate in situ, and the result is used to 
tunc 7r° kinematics in our simulation. We also used the measured ir° production 
uncertainty in our simulation. 46 

Another major background is the intrinsic background, namely v e (& e ) as beam 
contamination. Although they can be predicted by the beamline simulation, and 
are expected to be <0.5%, this is a critical background for the ^0.5% appearance 
oscillation search carried out by experiments, such as MiniBooNE. Again, in situ 
measurements largely help to reduce errors in the simulation. For example, the 
majority of v e (P e ) are from ^i-decay in the beamline but one can constrain their 
variations from the measured (y^) rate, where both v e and (y e and i> M ) are 
related through the n + (ir~) decay (for 7r + , 7r + — > ,u + — > 9 l _ l v e e + ). Another 

major source of v e (v e ) is kaon decay. MiniBooNE utilizes SciBooNE experiment 
data to constrain it. 47 SciBooNE is a tracker for the neutrino cross section mea- 
surement, located upstream of MiniBooNE, and their precise track measurement 
is sensitive to primary mesons in the beamline. More specifically, K-decay origin 
neutrinos are higher energy, and tend to make multiple tracks in the SciBooNE 
detector. This information provides the constraint on the errors on predictions of 
v e (v e ) from K-decay in MiniBooNE. 

After the evaluation of all backgrounds, MiniBooNE finds a signal-to- background 
ratio of roughly one to three, with expected oscillation parameters. 

4.5. MiniBooNE neutrino mode oscillation result 

For neutrino mode data analysis, 24 we use 6.46 x 10 20 protons on target (POT) data. 
After all cuts, an excess of v e candidate events in the "low-energy" region (200< 
£^ B (MeV)<475) was observed (Fig. 7). A total of 544 events are observed in this 
region, as compared to the predicted 409.8 ± 23.3(stat.)±38.3(syst.). Interestingly, 
this excess does not show the expected L/E energy dependence of a simple two 
massive neutrino oscillation model. Therefore, this excess might be new physics. 

4.6. MiniBooNE antineutrino mode oscillation result 

For the antineutrino mode analysis, 25 we use 5.66 x 10 20 POT data. Here, Mini- 
BooNE not only observed an excess in the low energy region, an excess in the 
"high-energy" region (475< £^ B (MeV)<1300) was also observed (Fig. 8). There- 
fore, in the "combined" region (200< £^ B (MeV)<1300), MiniBooNE observed 241 
v e candidate events as compared to the predicted 200.7 ± 15.5(stat.)±14.3(syst.). 
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Fig. 7. MiniBooNE neutrino mode u e appearance search result. The data points with errors are 
shown together with predicted backgrounds. 
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Fig. 8. MiniBooNE antincutrino mode i/ e appearance search result. The data points with errors 
are shown together with predicted backgrounds. 

Again, the i/SM does not predict this excess, which therefore has the potential to 
be new physics. 

5. Lorentz-violating neutrino oscillation analysis in MiniBooNE 

In this section we follow the procedure described in Sec. 3. 



5.1. The coordinate system 

First, we make the time distribution of neutrino events from a standard GPS time 
stamp. The analysis is based on the sidereal time distribution, but we also use 
the local solar time to check time-dependent systematics. The coordinate system is 
described in Fig. 3. The local coordinates of the BNB are specified by three angles, 
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local solar time (sec) sldeal time (sec) 



Fig. 9. Time distribution of the MiniBooNE low energy excess in neutrino mode. The top left 
figure shows the day and run number of each event. The bottom left plot is the scatter plot of each 
event with sidereal and local solar time; projections of the distribution onto each axes are show to 
the right. 

the co-latitude \ — 48.2°, the polar angle 9 = 89.8°, and the azimuthal angle 
<p = 180.0°. 

5.2. Time dependent systematics 

The signature of the Lorentz violation in this analysis is the sidereal time-dependent 
neutrino oscillations. However, any solar time distribution may show up in the 
sidereal time distribution, if data taking is not completely continuous. This is the 
case for MiniBooNE, as you see from top left plots of Figure 9 and 10. There were 
occasional shutdowns of the accelerator, and data taking was not continuous over 
any of the years. Therefore, we need to check time-dependent systematics, which 
include any day-night effects of the detector and the beam. For example, electronics 
may have a higher efficiency in cold nighttime than hot daytime, etc. The quick way 
to evaluate all of these effects together is to utilize the neutrino data itself. The high- 
statistics j^CCQE 35 data are used to evaluate day- night variations of the neutrino 
interaction rate. Figure 11, above, shows the variation of i^CCQE events in local 
solar time. 1 As you see, the j^CCQE sample shows day-night variations. These 
originate from the variation of the beam. The number of POT shows sinusoidal 
curve, namely maximum POT is at midnight and decreases towards the daytime: 
presumably the beamline is occasionally accessed in the daytime, and cumulative 
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Fig. 10. Time distribution of the MiniBooNE excess in antineutrino mode. The keys are same 
with Figure 9. 

human behavior over several years make a nice sinusoidal curve in POT, and in 
neutrino rate! After correcting POT, t^CCQE event shows flat (Fig. 11, bottom). 

Then, the question is how much the solar time variation of the POT actually 
shows up in the oscillation candidate sidereal variation. We found that a small 
variation persists in the f^CCQE sidereal distribution. We evaluated the impact of 
this on our analysis by correcting the POT variation event by event in u e candidate 
data. It turned out that, because of the low statistics of oscillation candidate events 
and the smearing effect from solar time distribution to sidereal time distribution, 
the correction only had a negligible effect. Thus, we decided to use uncorrected 
events in later analysis. 

Similar study is needed for the antineutrino mode, but the conclusion is likely 
to be the same due to lower statistics of D e sample. 

5.3. Unbinned likelihood fit 

To find best-fit (BF) parameters to describe MiniBooNE v e (v e ) candidate data with 
Lorentz violation, we employed an unbinned likelihood fit. The likelihood function 
A has following expression with two probability density functions (PDFs). 

A = Tn(^ + ^)x^e x p(-<^£) ,11) 
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Fig. 11. The top histogram shows the i/ M CCQE local solar time distribution, with the arbitrary 
normalized fit curve extracted from the POT local solar time distribution in the same time period 
(solid curve). The bottom plot shows the same events after correcting for variations in POT. The 
X 2 of this distribution with a flat hypothesis is 53.9/49 (29.3% compatibility). 



N, the number of observed candidate events 

(i s , the predicted number of signal events, the function of fitting parameters 

jib, the predicted number of background events, floating within la range 

F s , the PDF for the signal, the function of sidereal time and fitting parameters 

Tb, the PDF for the background, not the function of the sidereal time 

(Tfc, the lcr error on the predicted background 

fib, the central value of the predicted total background events 

This method is suitable because it has the highest statistical power for a low- 
statistics sample. The computation is performed to maximize this function. But in 
the reality, we maximize the log likelihood function. The maximum log likelihood 
(MLL) point provide the best fit (BF) parameter set. Then the constant surface of 
log likelihood function provide the errors. Neither the neutrino nor the antineutrino 
mode data allow us to extract errors if we fit all five parameters at once (Eq. 5), 
due to the high correlation of parameters. Therefore, we focus on three-parameter 
fit (Eq. 9) to discuss errors and limits. Since the five-parameter fit is quantitatively 
similar to the three-parameter fit, we will focus the discussion of these results on 
the three-parameter fit. 

6. Results 

This section describes the results of the fits. 
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Fig. 12. The fit results for the neutrino mode low energy region. The plot shows the curves 
corresponding to the flat solution (dotted), three- parameter fit (solid), and five-parameter fit (dash- 
dotted), together with binned data (solid marker). Here the fitted background is shown as a dashed 
line, however the best fit for the background is 1.00 (the default prediction). 



6.1. Fit result of neutrino mode data 

Figure 12 shows the neutrino mode low energy region fit results. 1 The solid and 
dash-dotted lines are the best fit curves from three- and five-parameter fit. The 
dotted line is the flat solution. Since the fit is dominated by the (C) efi parameter 
(sidereal time-independent amplitude) , both three- and five-parameter fit solutions 
have small time-dependent amplitudes and look like the flat solution. The details of 
the fit results were tested by using a fake data set. We constructed a fake data set 
with signal to evaluate errors of fit parameters of three-parameter fit ((C) eAI , (A s ) e ^n 
and (Ac)ep.)- Then we constructed a fake data set without signal, to evaluate the 
compatibility with flat solution over three parameter fit solution by A\ 2 method. 
It turns out data is compatible with flat solution over a 26.9%, and it concludes v e 
candidate data is consistent with flat. 



6.2. Fit result of antineutrino mode data 

Figure 13 shows the fit result for combined energy region in antineutrino mode, 
which is analogous with Figure 12. 1 For antineutrino mode, the combined region is 
more interesting due to lower statistics. The fit result is more curious in antineutrino 
mode, because now the (C) efl parameter no longer significantly deviates from zero, 
but the fit favors a nonzero solution for the (A s ) ef i and (*4 c ) eAI parameters. The 
fit solutions look more different from the flat distribution. We again constructed 
a fake data set to find the significance of this solution, and it turns out that the 
compatibility with the flat solution is now only 3.0%. Although this is interesting, 
however the significance is not high enough to claim discovery. 
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sidereal time (sec) 

Fig. 13. The fit results for the antincutrino mode combined region. Notations are the same as 
Fig. 12. The fitted background is shown as a dashed line, and the best fit for the background is 
0.97 (3% lower than the default prediction). 

6.3. Summary of the fit 

In this analysis, three parameters, (C) eM , (A s ) eu , and (A c )eu f° r neutrino and an- 
tineutrino mode are obtained by fits. Fits provide the BF values of above three 
parameters, as well as la errors and 2a limits. The expressions of these parameters 
are found in Eqs. 6, 7, and 8. From the 2a limits we obtain we estimate the limit of 
each SME coefficient by setting all but one of SME coefficient as nonzero. Table 1 
is the result. As you see, the limits of the SME coefficients from the MiniBooNE 
data are of the order of lO" 20 GcV (CPT-odd), and lO" 20 to 10" 19 (CPT-even). 
Similar analysis have been done for the LSND data. 49 However, these limits exclude 
any SME coefficients needed to explain the LSND data. Therefore, a simple picture 
using Lorentz violation to explain both LSND and MiniBooNE leaves some tension, 
and a mechanism to cancel the Lorentz- violating effect at high energy 17 ' 18 ' 19 ' 20 is 
needed. 

Conclusion 

Lorentz and CPT violation is a predicted signal at the Planck scale, and there 
are worldwide efforts to search for it. Neutrino oscillation is a natural interfer- 
ometer, and the sensitivity to Lorentz violation is comparable with precise opti- 
cal experiments. The MiniBooNE short-baseline neutrino oscillation experiment at 
Fermilab observed a v e (p e ) candidate event excess from a (v^) beam. These 
excesses are not understood by ^SM; therefore, they might be the first signal of 
new physics. Since Lorentz- violation-motivated phcnomenological neutrino oscilla- 
tion models, such as the Puma model, predict an oscillation signal for the Mini- 
BooNE, it is interesting to check the sidereal variation of MiniBooNE oscillation 
candidate data. The analysis found that neutrino mode data prefer a sidereal time 
independent solution, and the data is compatible with a flat distribution over 26.9%. 
The antineutrino mode prefers a sidereal time dependent solution, and the data is 
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Table 1. List of SME coefficient limits, derived from 2<x limits of fitting parameters, setting 
all but one of the SME coefficients to be zero. 



Coefficient c/i (u mode low energy region) c/j, {y mode combined region) 



Rc(cil) j or lm(<2i J ) J 


4.2xl0- 2u 


GeV 


2.6x10-™ 


GeV 


Re(a L ) x or lm(a L ) x 


6.0xltr 20 


GeV 


5.6xl0- 20 


GeV 


Re{ai,) Y or lra{ai,) Y 


5.0xl0- 20 


GeV 


5.9xl0- 20 


GeV 


Re(ai) z or lm(a L ) z 


5.6X10- 20 


GeV 


3.5xl0- 20 


GeV 


Rc(c L ) XY or lm(c L ) XY 










Rc(c L ) xz or Im(c L ) xz 


1.1x10" 


19 


6.2x10- 


20 


Rc(c L ) YZ or Im(c L ) YZ 


9.2x10" 


20 


6.5x10- 


20 


Ke(c L ) xx or lm(c L ) xx 










Rc{c L ) YY or lm(c L ) YY 










Rc(c L ) zz or Im(c L ) zz 


3.4x10" 


19 


1.3x10- 


19 


Rc(c L ) TT or lm(c L ) TT 


9.6x10" 


20 


3.6x10- 


20 


R c ( c l) TX or lra(c L ) TX 


8.4x10- 


20 


4.6x10- 


20 


Rc(c L ) TY or lm(c L ) TY 


6.9x10- 


20 


4.9x10- 


20 


Rc(c L ) TZ or lm(c L ) TZ 


7.8x10- 


20 


2.9x10- 


20 



compatible with a flat solution only 3.0%, making this solution very interesting; 
however, the statistical significance is not high enough to claim as evidence. Since 
the data set we used for this analysis is about ~50% of the total antineutrino mode 
data set, reanalysis including full data set may increase the significance. Finally, 
from the fits, we extract limits of each minimal SME coefficient. The results from 
MiniBooNE leave tension with LSND under the simple Lorentz violation motivated 
model. 
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